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A recent dynamical formulation at derivative level d'^g for fluid spacetime geometries {M, g, u), 
that employs the concept of evolution systems in first-order symmetric hyperbolic format, implies 
the existence in the Weyl curvature branch of a set of timelike characteristic 3-surfaces associated 
with propagation speed | | = | relative to fluid-comoving observers. We show it is the physical 
role of the constraint equations to prevent realisation of jump discontinuities in the derivatives of 
the related initial data so that Weyl curvature modes propagating along these 3-surfaces cannot be 
. . . activated. In addition we introduce a new, illustrative first-order symmetric hyperbolic evolution 

I system at derivative level d^g for baryotropic perfect fluid cosmological models that are invariant 

■ under the transformations of an Abelian G2 isometry group. 
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^ ! I. INTRODUCTION 

(N . 

^ ■ It is well known that the mathematical nature of the evolution system within the relativistic gravitational field 
CO ' equations is essentially hyperbolic, with domains of dependence and influence determined by the speed of light. 
However, this nature is not obvious when the dynamical equations are written out in standard form, employing either 



■ a metric approach , a Hamiltonian (ADM) representation [ 1 



o 
o 



^1; 



a 1 -I- 3 orthonormal frame (ONF) formulation , 
Consequently, over the years considerable effort has 



or a form obtained from a covariant (1 -I- 3)-decomposition 

gone into determining ways of making this hyperbolic nature clear. Indeed, evolution systems of partial differential 
\ equations in first-order symmetric hyperbolic (FOSH) format have proven to be a main theme of research for at least 
the last five years. This activity was particularly motivated by the desire to carry over to the numerical investigation 
of relativistic effects such as generation of gravitational radiation associated with the in-spiralling of black hole and 



I . neutron star binaries the methods and expertise gained in areas of computational physics with a longer tradition like, 
^JT)' e.g. , hydrodynamics, where evolution systems of FOSH format are commonplace (see, e.g., the recent reviews [ p4[ and 
• • [|13|| ) . 

One of the most promising methods for obtaining a FOSH representation for the relativistic gravitational field 
equations is to use an extended 1-f 3 orthonormal frame formulation that includes the once-contracted second Bianchi 
identities In terms of the highest-order partial derivatives of the spacetime metric g implicitly occurring, the 

■ dynamical equations here rank at derivative level d^g. The set of geometrically defined field variables contains the 
components of the physically significant Weyl curvature tensor, and so we refer to this as a 1 -I- 3 orthonormal 
frame curvature formulation of gravitational fields. It represents a completion of the physically and geometrically 
transparent 1 -f 3 covariant formalism for the cosmological case, in which all field variables are covariantly defined 
relative to a uniquely defined future-directed timelike reference congruence given by the unit 4-velocity field u of 
matter present As Friedrich first showed, when the matter source for a spacetime geometry (Al, g, u) is 

described phenomenologically as a perfect fluid, or one restricts to vacuum situations, upon introduction of a set 
of local coordinates and a specific choice to remove the gauge fixing freedom it is indeed possible to find linear 
combinations of the field variables and their dynamical equations that lead to an evolution system of (autonomous) 
partial differential equations in FOSH format, as desired ||ll| , |l^ . It is puzzling, then, that when this is done, as well as 
the expected sets of characteristic 3-surfaces associated with propagation speeds relative to u of | u | =0 (Coulomb-like 
gravitational and fluid rotational modes), \ v\ — Cg (sound wave modes), and | w | = 1 (transverse gravitational wave 
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modes), the equations on the face of it indicate that the semi- longitudinal Weyl curvature characteristic eigenfields, 
identified long ago by Szekeres , appear related to an additional set of timelike characteristic 3-surfaces associated 
with propagation speed relative to u of one- half the speed of light, \ v\ = ^. 

Friedrich assumed that these modes could be ignored, because, at least in the vacuum case, they had no invariant 
meaning Two of us (HvE and GFRE), on the other hand, argued in a recent paper that in the fluid context 

that is relevant in relativistic cosmology, this is not the case, essentially because (i) the timelike reference congruence 
corresponding to u is invariantly defined by the presence of matter, and (ii) the hypothetical Weyl curvature modes 
cannot all be transformed to zero by using up the remaining gauge fixing freedom ||^. Thus, these might possibly 
represent propagation modes that are physically meaningful. In private communications, Friedrich responded that 
the form of the d^g-oidei FOSH evolution systems discussed in Refs. and Q is not unique; one could obtain 
other reductions of the relativistic gravitational field equations that would lead to additional timclike characteristic 
3-surfaces associated with propagation speeds different from \ v\ = i. And, hence, these modes cannot be physical 
(see also Ref. [Q). However, that argument is not very satisfactory. In our view, it compounds the problem, rather 
than solving it, by suggesting (if we take for face value the FOSH formalism and its interpretation) that still other 
(subluminal) propagation speeds could be associated with the transport of physical quantities, as well as \ v\ = ^. 
The implication is that any single FOSH evolution system cannot be taken seriously without considering the set of 
all possible FOSH evolution systems, which is difficult to determine. 

What had not been fully taken into account in these discussions, even though Sec. Ill D 2 of Ref. [|j made a brief 
reference to this aspect, is the physical role the constraint equations play in the selection of appropriate classes of 
initial data sets. Specifically, what kind of physically relevant jump discontinuities they do allow to be present in the 
derivatives of these data sets ||2^. Jump discontinuities in the derivatives of the initial data are of major dynamical 
importance as they can be thought of as representing the triggering or terminating events of physical generation 
processes.^ In the present paper we show that an analysis of the constraint equations leads to a resolution of the 
problem: while the additional timelike characteristic 3-surfaces do indeed occur in the d^g-order FOSH evolution 
system, the constraint equations are of such a form as to prevent the corresponding semi-longitudinal Weyl curvature 
modes being activated. Because jump discontinuities cannot exist in the values of the constraint equations themselves 
(they have to be zero everywhere), jump discontinuities in the derivatives of the initial data can only be propagated 
along the characteristic 3-surfaces associated with propagation speeds relative to u of 1 1; | =0, | w | = Cg, and \ v\ — 1, 
but not along those associated with \ v\ = i. Hence, one cannot physically send arbitrary information along the latter 
set of timelike characteristic 3-surfaces. 

The problem of additional characteristic 3-surfaces and how the constraint equations select appropriate initial data 
sets can be highlighted in a nice and transparent fashion by a suggestive example that derives from a set of linearised 
relativistic gravitational fi eld equations given in a paper by Kind, Ehlers and Schmidt p^ . We will briefly discuss 



this example in Appendix A 1 



A possible complementary explanation to the above viewpoint arises from the following chain of arguments: 

(i) For vacuum as well as perfect fluid spacetime geometries, there exist c^^g-order dynamical formulations of the 
relativistic gravitational field equations (with harmonic coordinate gauge fixing) in which characteristic 3-surfaces 
only correspond to the local light cones, the local sound cones, or are generated by the fluid flow lines. (See, e.g., 
Refs. [H and 11.) 

(ii) There are dynamical formulations of order d^g or d^g that have additional (timelike or spacelike) characteristic 
3-surfaces. Along these characteristic 3-surfaces the evolution system will propagate jump discontinuities if we put 
them into the related initial data {ignoring the constraint equations). 

(iii) Different dynamical formulations with different sets of characteristic 3-surfaces all describe for identical initial 
data the same spacetime geometry. 

(iv) Because of the existence of a dynamical formulation according to (i), jump discontinuities as described in (ii) 
cannot really exists (and so be physically meaningful). Hence, the only logical possibility is that in (ii) the constraint 
equations prevent the existence of such jump discontinuities. 

The present paper shows that the semi-longitudinal Weyl curvature modes that potentially exist cannot in fact carry 
"gravitational news" . Of course, this is expected, so this analysis simply confirms what everyone has believed all along. 
However, we believe it is nevertheless useful, both in terms of showing important relations between the constraint 
equations and the set of characteristic 3-surfaces that can occur in general FOSH evolution systems (and have not 
been clearly demonstrated in the literature), and because this does indeed resolve the issue at hand in the specific 
important case of relativistic gravitation, where the variables presented here, and hence the associated dynamical 
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equations, have much to recommend them. Without the analysis given in this paper, the apparent (| w | = ^)-modes 
remain an annoying and unresolved problem. 

The outline of the paper is as follows. In Sec. |l| we briefly discuss the mathematical concepts relevant to our analysis: 
we review the central ideas behind FOSH evolution systems, we address the issue of gauge fixing freedom in a 1 + 3 
orthonormal frame formulation of the rclativistic gravitational field equations, we list the set of constraint equations we 
need to consider, and we introduce a (l + l + 2)-decomposition of all geometrically defined field variables. Then, in Sec. 



[II, we address from a purely local viewpoint the question of how the constraint equations select appropriate initial 
data sets for the d^g-oider FOSH evolution system introduced in Ref. |9|. To the best of our knowledge, we present 



in Sec. IV for the first time a fully gauge-fixed autonomous FOSH evolution system for a special class of spatially 
inhomogeneous perfect fluid cosmological models on the basis of a 1 + 3 orthonormal frame formulation that is only of 
order d^g in the degrees of freedom of the gravitational fieldj^ We refer to this as a 1 + 3 orthonormal frame connection 
formulation of gravitational fields. Additionally, we derive for these cosmological models, for which there exists an 
Abelian G2 isometry group, the transport equations that describe how physically relevant jump discontinuities in 
the derivatives of the initial data are propagated along the bicharacteristic rays of the set ting . Our conclusions are 
contained in Sec. Finally, besides the above-mentioned suggestive example in Appendix Al , we give in Appendix 



A 2 in explicit form those linear combinations of the components of the extended 1-1-3 orthonormal frame constraint 
equations that prove suitable for our analysis. 

We will use the same conventions, units and notations as introduced in Appendix A 1 of Ref. 



II. MATHEMATICAL PRELIMINARIES 

A. First-order symmetric hyperbolic evolution systems 

We consider evolution systems for a collection of k real- valued field variables — u"^{x^) that are composed of a 
set of k quasi-linear partial differential equations of first order given by 

M^^^'(x^^^'^)5^^^B = , A, B, C = 1, . . . ,k ; (2.1) 

the field variables are functions of a set of local spacetime coordinates { x'' }. Evolution systems of this form are 
called symmetric if the real-valued k x k coefficient matrices entering the principle part satisfy M^^ = M^'^^^'^; 
moreover, they are called hyperbolic if the contraction M^^ with the coordinate components of an arbitrary past- 
directed timelike 1-form Ua yields a positive-definite matrix. We remark that (i) cases with M^^ ^ = M"^^ ^{x") are 
referred to as semi-linear , and (ii) cases with M^^ ^ = M^^^(it'^) and — N ^ju'^) are referred to as autonomous. 
In general, it proves convenient to consider a (1 -I- 3)-decomposition of Eq. ( ^.l| ) in the format 

M'^^°{t,x\u^)^tUB + M'^^'{t,x^,u^)^^UB-^N'^{t,x\u^) . 

The concept of FOSH evolution systems was first introduced by Friedrichs ||l^; a standard reference, broughtly 
discussing its usefulness to applications in mathematical physics and also presenting proofs for existence and uniqueness 
of solutions, is the book by Courant and Hilbert The characteristic condition 

= Q:=dct[M'^^'^V^(/)] (2.2) 

determines the coordinate components of the past-directed nor mals Va4' oi the set of characteristic 3-surfaces 



C:{(j){xJ^^ const} associated with the FOSH evolution system (|j). With M^^'' = M^-^^^^, hyperbolicity of 
Eq. ( |2.1| ) thus also corresponds to all individual roots ("eigenvalues") v of Eq. (2.2) being real-valued. Every 
individual v then defines a pair of so-called left and right nullifying vectors, l^ and r^, by 

0^:1 A (Af-^^ ^ V^c/.) , (Af^^ ^ V^c/.) rs ; (2.3) 

the linearly independent sets {l^} or { 7"^ } form a basis of the fc-dimensional space of field variables u"^. 



■^The equations of motion for the matter sources will necessarily derive from the twice-contracted second Bianchi identities 
and so are typically of order d^g. 
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A ccording to the theory discussed in Ch. VI. 4. 2 of Courant and Hilbert ||], FOSH evolution systems of the format 
( ^.l|) have the power to describe the physical transport along so-called bicharacteristic rays of jump discontinuities 
that exist in the outward first derivatives across a characteristic 3-surface C:{(j){x'^) = const} of the field variables 
u^; the tangential first derivatives of the as well as the themselves are assumed to be continuous across 
C:{(/)(a:'') — const}. As is standard, we will use the notation 

[/] := hm / - hm / = /+ - 

(p — >C-^- q> — >C — 

to symbolise a jump discontinuity (of finite magnitu de) acros s C:\ (j}{x^^) — const} in the value of a given variable /. 



Under the stated assumptions, it follows from Eqs. {2A) and ( |2.3| ) that 



= {M^B ^ V^0) [ d^UB ] ^ [ d^u^ ] = [d^u] r^ , (2.4) 

i.e., the jump discontinuity [S^m"^] must be proportional to a right nullifying vector r^ . The real- valued scalar of 
proportionality, denoted by [S^m], is assumed to have continuous first derivatives. Then, according to Chs. VI. 4. 2 
and VI. 4. 9 of Ref. |^, for linear, semi- linear and quasi-linear FOSH evolution systems (|]^) the transport equation for 
[d^u] along bicharacteristic rays within the characteristic 3-surfaces C:{(l){x'^) — const} takes the effective form 

= {lAM^^'^^rB) d^[d^u] + ( (Z^Af^^^) d,,rB - {IaN^'^tb) ) [B^u] . (2.5) 

Note, in particular, the involutive character of this relation; if [drj,u\ is non-zero at one point along a bicharacteristic 
ray, it will be non-zero everywhere along this ray, and vice versa. Note also that the present treatment of jump 
discontinuities breaks down when the C:{0(x^) = const} within a given family start to intersect and so pro mpt the 



formation of so-called "shocks". Shock formation, however, cannot arise when the principal part of Eq. (2.1) is 
semi-linear. It is a special feature of the relativistic gravitational field equations that related FOSH evolution systems 
do have, in the branches that evolve the degrees of freedom in the gravitational field itself, principal parts which are 
effectively semi-linear.^ More precisely, the coefficient matrices M'^^^ in these branches depend only on those field 
variables which form a background for gravitational dynamics in the sense of Geroch p5[ . 

B. Choice of gauge source functions and local coordinates 

As Friedrich emphasised in Sec. 5.2 of Ref. |ll[, there exists within a 1 -I- 3 orthonormal frame representation of the 
relativistic gravitational field equations a set of ten so-called gauge source functions, G := { T°, T", T"o, T"fj }, that 
can be arbitrarily prescribed in any dynamical consideration (and are thus assumed to be "known"). These relate to 
(i) the arbitrary choice of a future- directed reference "time flow vector field" T,^ which, in terms of the 1 -f 3 ONF 
basis { eo, e^}, is expressed by 

T :=T°eo + T"e„ , T" > , (2.6) 

and (ii) the propagation of the 1 -f 3 ONF basis { eo, } along T, described by 

VTeo:=r"oe„, Vxe^ := Gq + T'^^ . (2.7) 

Parallel transport of { ep, } along T, for example, thus corresponds to setting — T^a — T^p. Upon introduction 
of a dimensionless local time coordinate t and dimensionless local spatial coordinates { x' } that comove with T, the 
gauge conditions related to a 1 -f 3 orthonormal frame representation are made explicit by | [ll| ] 

eo'^j;^{M^'S%-T'^e^n, r"„o = ^ (T"a - r% T'^) , po ^ ^ {T" p - p.T^) ; (2.8) 

we keep the inverse unit of [length], Mq^, as a coefficient for reasons of physical dimensions. 

The fluid-comoving, Lagrangean perspective adopted in the discussion of Ref. ^ by identifying the timelike reference 
congruence with the fluid 4- velocity field, eg = u, is now obtained by fixing three of the four dimensionless coordinate 
gauge source functions according to = 0, resulting in an alignment T || eg (= u). This leads to 



^See, e.g., Eqs. (3.21) - (3.23) in Ref. g, or Eqs. ( [l.l7| ) - ( [l.20| ) in the Abelian G2 example given in Sec. |V 
*The reference vector field T need not necessarily be timelike. 



below. 
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eo'' = M-'S%, =T°r^o = T°t(„ , T"^ = r° T^^o = T° e"^^ fl^ , (2.9) 

where M := T^Mq. Consequently, the three frame gauge source functions T^a become proportional to the components 
of the fluid acceleration u", while the three frame gauge source functions T"^ become proportional to the components 
of the rotation rate fi" at which the spatial frame { Go. } fails to be Fermi-propagated along u. For cosmological models 
(tW, g, u) with perfect fluid matter sources, Ref. g] introduced proper time along u by setti ng T " — 1 ^ M = Mq, 
and derived the evolution equation for along u from the commutators on the basis of Eqs. ( ^.12 ) and (2.21) below. 



Additionally, Ref. ||9[] set fJ" = 0. The latter choice, however, is by no means compulsory, and other choices may prove 
equally convenient (given that Vl" is assumed to be "known" ) . 

In order to obtain from the extended 1 + 3 orthonormal frame relations proper partial differential equations such 



that the theory underlying Subsec. II A applied, Ref. expressed the coordinate components eo'' := eo(x^) and 
Bq,^ := ea{x^) of the 1 + 3 ONF basis { Gq, e^, } in terms of the comoving local coordinate basis { 9t, di } by 

Go Af-i dt , G„ {M, dt + d,) ; (2.10) 

M — M[t,x^) denotes the threading lapse function and Midx^ — Mi(t, x^) dx'' the dimensionless threading shift 
1-form. The inverse of the threading metric is h^^ := 5°^^ Bq* e^^ . See, e.g., Ref. and references therein. 

C. Matter model 

The matter sources in the present discussion (as well as in Ref. |^ ) are assumed to be modelled phenomenologically 
as a perfect fluid such that, with respect to fluid-comoving observers, 

= (7"(u) = 7r„/3(u) , (2.11) 

i.e., the energy current density and the anisotropic pressure both vanish. Additionally, a baryotropic equation of state 
is assumed, 

p = p{p), (2.12) 
relating the isotropic pressure p(u) to the total energy density /i(u). 

D. Constraint equations 

The following relations in the set obtained from an extended 1+3 orthonormal frame representation of the relativistic 
gravitational fleld equations do not contain any frame derivatives with respect to gq. Hence, it is commonplace to 
refer to these relations as "constraint equations".^ These are |19 1^ 



= (Ci)" := (G;3-3a^)(a"^)-|5"'5G/3(e)-n"^c^^ + e"''^[(G;3 + 2u0-a^)(u;^)-n05(T*^] (2.13) 

= (C2) := (G,-?i„-2a„)(w") (2.14) 

= (C3)"^ := ((ST<"G^ + 2u<" + a<")(a;'^>)-in\a"'^ + 3n<"^f7^>^ + F"'' 

- e^*<" [ (g^ - a^) (af^h) + u;5 ] (2.15) 

= (Cj)" := (6^-2a0)(n"'3) + |ecj" + 2(T";3Cc;^ + e"^^[6^(a^)-2c<j^f]^ ] (2.16) 

= (Cg)"'' := *S'"'^ + ie(T"''-cr<"^cr'5>'^-w<"w'3> ^2w<"f7'5) (2.17) 



^Even though this terminology is problemat ic in the generic case when eo = u has non-zero vorticity, <-o"{u) ^ 0, and local 
coordinates are introduced according to Eq. (p.l[)|) above. 
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where 



= (Cg) := *i? + I - {(To^po"''^) + 2 {uo,uj°') - 4 f]") - 2 /I - 2 A 

= (Cs)" := - 3 ap) {H'^'^) + {y. + p) uj"^ + i lo p E^" + e"'^'' [ aps E% ~ nps H% ] 
= (Cpf)" := clS''^ep{pi) + {pi+p)u'' , 

*Sap := e(a(a^)) + b(^ap) - (a (e|-y| - 2a|^|) (n^)^) 



*i?:=2(2e„-3aa)(a")-i6"„ 



(2.18) 
(2.19) 
(2.20) 
(2.21) 

(2.22) 
(2.23) 



2n, 



a7 



(2.24) 



c^(/i) := dp{y)/dy defines t he ise ntr opic speed of sound with < < 1, and angle brackets denote the symmetric 



tracefree part Q While Eqs. (^.13|) - (|2.18|) derive from the Ricci identities, the Jacobi identities and the Einstein field 
equations and so are of order d'^g, Eqs. (2.19) - (2.21) derive from the o nce- and twice-contracted second Bianchi 
identities and so are of order d'^g. The divergence constraint equation ( 2.13 ) for the fluid rate of shear is often referred 
to as the "momentum constraint". When = cj"(u), such that th e flu id 4- velocity field u constitutes the normals 
to a family of spacelike 3-surfaces S:{t = const}, Eqs. ( ^.17 ) and ( 2.1g ) correspond to the symmetric tracefree and 
trace parts of the once-contracted Gaufi embedding equation. In this case, one also speaks of (Cg) as the generalised 
Friedmann equation, alias the "Hamiltonian constraint" or the "energy constraint" . 



E. (1 + 1 + 2)— decomposition 



In the present discussion it proves very helpful to consider a (1 + 1 + 2) -decomposition of all geometrically defined 
field variables and their dynamical relations. In order to do so, we arbitrarily pick the frame basis field ei as a second, 
spacelike, reference direction, in addition to eg = u as a timelike one; any other spatial direction, however, would be 
equally acceptable. Hence, in a small isotropic neighbourhood U in the local rest 3-space of an arbitrary event V, we 
establish the convention of regarding those spatial frame components of geometrical objects which contain the index 
"1" as (semi-) longitudinal with respect to ei, while regarding those which exclude the index "1" as transverse with 
respect to ei. Likewise, in U, ei shall constitute the outward frame derivative while 62 and shall be tangential 
frame derivatives when we turn in Sec. Ill to the issue of jump discontinuities across spherical spacelike 2-surfaces 



J^:{t = const, 4>{x^) — const}. For the frame components of spatial rank-2 symmetric tracefree tensors Oap — 
with squared magnitude :— ^{oapa^f^) > 0, we define a new set of frame variables by 



a+ := 2 (a22 + 033 j 
1 



a 



V3 



0.23 



: an a- := 



02 



7= (a22 - 033) 



2V3 
1 



V3 



031 



°3 •- 73 ""l^ - 



SO that 



3 ("4 



(2.25) 



(2.26) 



In particular, in the present discussion we have Oap € { (Tap, Eap, HapV. We remark that these definitions are now 
adapted to the conventions of the book edited by Wainwright and Ellis , implying they differ by a factor of -i from 
those used in Ref. ||^. In analogy to Eq. ( ^.25 ), we perform a (1 -I- 1 -|- 2)-decomposition of the spatial commutation 
functions Uap by defining 



The numbering of the constraint equations we employ is based on the conventions established in 1 + 3 covariant treatments 
of relativistic cosmological models {M, g, u) (cf. Ref. Q). 
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n := nil + 7122 + "33 := - "ii + 5 ("-22 + 'T'Sa) JT-- := (n-22 - "-33) 

(2.27) 

"-x:=^"23 ri2:=^n3i ^^3 ^^l2 • 

The squared magnitude is then given by 

i(n„0n"'3) - i (n^ + 2nl) + 3 {nl + +71^ + n^) . (2.28) 

Note that only (n — 2n+), 7i_ and transform as tensor components under rotations of the spatial frame {e^} 
about the reference ei-axis. 



By employing these conventions and d efinit ion s, we have listed in Appendix A 2 certain linear combinations of the 



components of the constraint equations ( 2.13 ) - ( 2.21 ) that will be needed in the following sections 



INSERT Table | on "Conventions for (1 + 1 + 2)-dccomposition" HERE. 



III. PHYSICAL EFFECT OF CONSTRAINT EQUATIONS ON OUTWARD FIRST DERIVATIVES 

Generic cosmological models ( A^, g, u) with a perfect fluid matter source have fluid 4- velocity flelds u with non- 
zero vorticity, lo°'{\i) ^ 0. This property makes it impossible to determine a fluid-comoving spacelike 3-surface 



S:{t = const} everywhere orthogonal to u on which initial data satisfying the constraint equations of Subsec. II D 
could be specified. In this case, the discussion of a wcU-posed Cauchy initial value problem requires that the setting 
of the data as well as the solution of the constraint equations be instead performed on a non-comoving spacelike 
3-surface S:{t = const}, before the data is evolved along u.|^ All these complications disappear when — uj"{u), and 
well-defined spacelike 3-surface S:{t = const} everywhere orthogonal to u do exist. 

For our purposes, however, it is sufficient to investigate the physical effect of the constraint equations on the 
selection of appropriate initial data sets from a purely local viewpoint, i.e., only in a small isotropic neighbourhood U 
in the local rest 3-space of an arbitrary event V. Due to the local Minkowskian structure of all relativistic spacetime 
manifolds g ) , one conventionally determines (and analyses their physical properties) the set of characteristic 
cones C:{4>{x'^) = const} for a given FOSH evolution system only within the small isotropic spacetime neighbourhood 
{— e < t < e} X U oi V. In line with this, we will consider in the following spherical spacelike 2-surfaces J^:{t = 
const, (/)(2:^) = const} in U across which we assume to exist (i) jump discontinuities in the outward first frame 
derivatives of certain geometrically defined field variables W^, and (ii) continuity of the tangential first frame derivatives 
of the and the themselves. As the constraint equations have to be satisfied ever ywh ere, i t is clear that across 



J' we have = [(C)"' ' Jt.^^const for the value of any component in the set of Eqs. ( 2.13 ) - ( 2.21 ). 

Motivated by the prospect of grounding the discussion on wave-like phenomena described by the relativistic 
gravitational field equations on the deviation equation for a set of test particles, the dynamical considerations on 
the d'^g-oideT FOSH evolution system in Ref. ||] focused on the set of Weyl curvature characteristic eigenfields 
{ E+, H+, {E3 TH2), {E2 ± -ffa), {E- =F Hy_), {E^ ± i?-) }■ There it was correctly argued that the deviation equation 
monitors the physical effects on the state of motion of a set of test particles of both gradual as well as sudden changes 
in the values of these fields (see, e.g., Refs. and |26|) . What was overlooked in this work, however, is the fact 
that, on the basis of the theory underlying Subsec. p^I A| , the d^g-ovder FOSH evolution system presented in Ref. 
1^ can at best describe the physical transport along bicharacteristic rays of jump discontinuities in the (outward) 
first derivatives of these fields rather than these flelds themselves, given the constraint equations do not impose any 
additional restrictions. It is our aim to supplement the discussion of Ref. |^ by such a consideration in the present 
section. 

A. Jump discontinuities at derivative level d^g 

Considering in U a spherical spacelike 2-surface J^:{t — const, (l){x^) = const}, and assuming that across J' all 
geometrically defined field variables as well as their tangential first frame derivatives are continuous, we find that 



^It is currently unknown whether the d^g-order FOSH evolution systems with perfect fluid matter sources in Refs. and 
M can be generalised to a non-comoving perspective. 
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the Weyl curvature divergence equations (All) - (A14) amongst the constraint equations lead to the following set of 
jump conditions: 



[ei(^^+)]t,0=const 
[ei(£'3 T^^2)]t^0^const 

[ei(ii;2±i?3)],^^_t 

[ei(£;x±i/-)]*,^.eonst 







unconstrained 
unconstrained . 



(3.1) 
(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 



The implications are three-fold: 



(i) Jump discontinuities in the outward first frame derivative of the Coulomb-like Weyl curvature characteristic 
eigenfield originate from jump discontinuities in the outward first frame d eriva tive of the matter total energy 
density /i and are physically allowed, if the momentum conservation equation (A15) is satisfied on both sides of J 
with = [(CpF)i]t.</)=const- In more detail: assuming that all of p, Cg and u\ are continuous across J ^ Eq. ( A15[ ) yields 



= [(CpF)l]t,0=const = <?s [ei(Ai)] 



const 



const 



(3.7) 



Hence, if p = Cs = and ui = on both sides of J", the values of [ei(//)]t_0=const and [/Lt]t, ,/,=const remain 
unconstrained. Phenomena of the presently described kind occur, for example, across the surfaces of static, spherically 
symmetrical perfect fluid stars with equation of state (2.12) (see, e.g., Ref. [^). In a d^g-oxAex formulation, it follows 
that real- valued initial data for /i (and so for Ej^) is required to be of differentiability class C^{U) with respect to 
the zeroth-order derivative level of g. The equations of Ref. |^ show that [ei(i?+)]t,0=const ^ propagates with 
characteristic velocity v — relative to u. Note that in the vacuum subcase ^ = [ei (i?+)]t, ^^const = 0. Jump 
discontinuities in ei(_ff+) are not physically allowed, and so, in a 9'^g-order formulation, real- valued initial data for 
is required to be of differentiability class C^iU) with respect to the zeroth-order derivative level of g. 

(ii) Jump discontinuities in the outward first frame derivatives of the semi-longitudinal Weyl curvature characteristic 
eigenfields {E^^H2) and {E2±H^), that by restricting to the net d'^g-order FOSH evolution system in Ref. ^ {without 
accounting for the constraint equations) are theoretically associated with characteristic velocities w = ± ^ relative to 
u, are not physically allowed. Hence, in a 9^.g-order formulation, real-valued initial data for (£'3 =p H2) and {E2 ± H^) 
needs to be of differentiability class C^iU) (rather than C^iU)) with respect to the zeroth-order derivative level of g. 

In short, not all initial data can be given freely. 

It can b e easily inferred from the propagation equations along u for the Weyl curvature divergence equations ( 2.19| ) 
and ( ^.20 ), first published for an irrotational pressure- free fluid matter source in Ref. |18|, and presented for a general 
perfect fluid matter source in Refs. Q and that the characteristic velocities relat ive t o u fo r the components 
{Ci)2T{Cz)s, and ((74)3 ±((75)2 are w = ± ^ too. Hence, comparing this result with Eqs. ( A13| ) and (A14) in Appendix 
and Eqs. (3.3) and (3.4) above, it becomes clear that the Weyl curvature divergence equations propagate relative 



to u at precisely the speed that is required to ensure that jump discontinuities in ei(iJ3 ^ H2) and ei{E2 ± H^i) 
will always remain physically disallowed at any instant throughout the dynamical evolution of a cosmological model 
(Al, g, u). It should be emphasised at this st age that t his property is completely independent of the presence of 
matter. That is, of course the jump conditions ( |3.3| ) and ( |3.4| ) apply equally to vacuum spacetime configurations. 

(iii) Jump discontinuities in the outward first frame derivatives of the transverse Weyl curvature characteristic 
eigenfields [E- =F iJx) and (i?x ± H-) are physically allowed. Clearly, this situation reflects the freedom of specifying 
four arbitrary (non- analytic) real- valued functions Ig^g :— {ai(a;'), 02(2;*), 03(0;*), 04(0;*)} of differentiability class 
C^{U) with respect to the zeroth-order derivative level of g as the initial data for the dynamical degrees of freedom 
associated with the gravitational field itself. 



B. Jump discontinuities at derivative level d g 



To be able to argue in terms of physical effects described by the deviation equation for a set of test particles, we 
have to turn our attention directly to the set of Weyl curvature characteristic eigenfields and possible discontinuous 
changes in their values. Such changes are driven by the dynamics of the underlying connection fields and their 
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derivatives. Hence, to facilitate the interpretation of generic gravitational dynamics, it would be desirable to have 
available a d^g-oidei FOSH evolution system derived from a 1 + 3 orthonormal frame connection formulation of 
gravitational fields (see Refs. ||l9| and [|o| for reviews of the latter). Unfortunately, to date, such a formulation has 
not been accomplished for the generic case. The exception is the c^^g-order FOSH evolution system for perfect fluid 
cosmological models {M, g, u) with an Abelian G2 isometry group we will present in Sec. |^ below. In the absence 
of such a generally applicable dynamical formulation, we return to our local viewpoint and investigate how, in a small 
isotropic neighbourhood U in the local rest 3-space of an arbitrary event V, the constraint equations at derivative 
level d^g restrict the occurrence of jump discontinuities in the values of the Weyl curvature characteristic eigenfields 



themselves. To this end, we now focus on E qs. (^-13^ - (p-18 ), and certain linear combinations of the components 
thereof provided by Eqs. ( [A5| ) - ( [A10| ) and ( Aisf ) - (|A23[ ). Again, we consider in U a spherical spacelike 2-surface 



J:{t = const, — const}, and assume that across J' all geometrically defined field variables 

tangential first frame derivatives are continuous. This leads to: 

(i) w = longitudinal Weyl curvature characteristic eigenfields: 

= — i ei(ai) + tangential frame derivatives/algebraic terms (3.8) 
H+ = i ei((jji) + tangential frame derivatives/algebraic terms , (3.9) 



from Eqs. (AI8) an d (AT9| ). Across J, the generalised GauB-Friedmann equation (AlC) and the fluid vorticity 



divergence equation (A9), respectively, then impose the restrictions 

[-E'+]t,0=const = ^ 5 [^l("l)]t,0=const = Nt,</,=const (3-10) 

[^+]t,0=co„.t = \ [ei('^l)]t,0=const = • (3-11) 

That is, jump discontinuities in the values of the Coulomb-like Weyl curvature characteristic eigenfield i?+ originate 
from jump discontinuities i n the values of the matter total energy density /i and are physically all owe d, if the momen- 



tum conservation equation (|A15|) is satisfied on both sides of with = [(CpF)i]t,,/,=const [cf. Eq. (|3.7|)]. Ina9^g-order 
formulation, real- valued initial data for (and so for E+) is thus required to be of differentiability class C^iU) with 
respect to the zeroth-order derivative level of g. Note that in the vacuum subcase /i = [i?+]t, ,/,=const = 0. Jump 
discontinuities in the values of iJ+ are not physically allowed, and so, in a 9^g-order formulation, real-valued initial 
data for is required to be of differentiability class C^{U) with respect to the zeroth-order derivative level of g. 

(ii) w = ± i semi-longitudinal Weyl curvature characteristic eigenfields: 

{E3 =F H2) = T ^ ©1 (0-3 =F ?i2 " 073 =F ^2) + tangential frame derivatives/algebraic terms (3-12) 
{E2 ± H3) ~ T ^ ei((T2 ± 713 -I- a;2 T 13) + tangential frame derivatives/algebraic terms , (3.13) 



from Eqs. ( A2C ) and ( A21 ). In this case, we find that across J the fluid shear divergence/ Jacobi constraint equations 
(A7) and (|A8|) impose the restrictions 

[{E3 T ^2)]t^^^eonst = T 5 T n2 - ^ ^3 T ^ a2)]t>=co„«t = (3.14) 

[(£^2 ± ff3)]t,^=eonst = T ^ [ei(a2 ± 71^ + ^ ^2 T ^ a3)]t,0=co„st = . (3.15) 

That is, jump discontinuities in the values of the semi-longitudinal Weyl curvature characteristic eigenfields {E3 ^ H2) 
and {E2 ± H3) are not physically allowed. In a 9^5-order formulation real-valued initial data for (£^3 =F H2) and 
{E2 ± H3) is thus required to be of differentiability class C'^iU) (rather than C^{U)) with respect to the zeroth-order 
derivative level of g. 

(iii) V = ±1 transverse Weyl curvature characteristic eigenfields: 

(£'_ =F Hy.) ~ =F ei(fT_ Tn-x) + tangential frame derivatives/algebraic terms (3.16) 
(£'x ± — =F ei(fTx ± n_) + tangential frame derivatives/algebraic terms , (3.17) 



from Eqs. (A22) and (A23). In this case, we find that across the constraint equations impose no restrictions so 
that 

[{E- T £^x)]t,^=eonst = T [ei((T- T '^x)]t,^^const = uuconstraincd (3.18) 
[{Ex ± i?-)]t_^=eonst = =F [ei((Tx ± )] t,^=eonst = Unconstrained . (3.19) 
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That is, jump discontinuities in the values of the transverse Weyl curvature characteristic eigenfields =F ff x ) 
and {Ex ± H-) are physically allowed. Hence, they can transport arbitrary non-zero jump discontinuities [ei((T_ ^ 
rix )]t,0=coiist and [ei(crx ± "--)]*, 0=const of finite magnitude. Again, this reflects the freedom of specifying four arbitrary 
(non-analytic) real-valued functions Ig2g := {Ai(x'), ^2(0;*), A3(x'), ^4(3:;*)} of differentiability class C^iU) with 
respect to the zeroth-order derivative level of g as the initial data for the dynamical degrees of freedom associated 
with the gravitational field itself. 



For com pleteness, we now also briefly discuss the effect of the constraint equations ( ^.13 ) - ( 2.16| ), respectively, 



Eqs. (A5) - (|A9D, on the outward first frame derivatives of the characteristic eigenfields associated with the fluid 
kinematical branch of the d^g-ordei FOSH evolution system in Ref. ^ [cf. Eqs. (3.27) and (3.28)]. The derivative 
level is hence one below the Weyl curvature case. For the v — Q fluid kinematical characteristic eigenfields we find the 
jump conditions: 



(iv) w = fluid kinematical characteristic eigenfields: 

[ei(ie + a+)]t,^=const = (3.20) 
[ei(wi)]t,,/,=const = (3-21) 
[61(0-3 -f ^ W3)]i,0=const = unconstrained (3.22) 

[ei(cr2 - ^'^2)]t,0=const = unconstrained . (3.23) 

The last two conditions imply the existence of two generically non-zero fluid rotational modes that were identified 
before by Ehlers et al in an analysis of linearised perturbations of arbitrary background dust spacetimes [Q. Cor- 
responding real- valued initial data for these modes is required to be of differentiability class C^{U) with respect to 
the zeroth-order derivative level of g. Finally, for the different parts of the w = ± fiuid kinematical characteristic 
eigenfields we find the jump conditions: 



(v) V — ±Cs fluid kinematical characteristic eigenfields: 

[ei(i e - 2o-+)]t_0=const = unconstrained 

[ei(0'3 — ^3)]t,4>=const — ± [®l("73 "2 + '^2)]t,0=const 
[61(0-2 + ^ W2)]t^0=const = ± [ei(^ 03 - "3)] 

[ei(ui)]t^0=const = unconstrained 
[ei(u2)]t,0=const = unconstrained 



[ei(u3)]t,0=c 



— unconstrained . 



= unconstrained 
=const = unconstrained 



(3.24) 
(3.25) 
(3.26) 
(3.27) 
(3.28) 
(3.29) 



Again, corresponding real- valued initial data for these parts is required to be of differentiabil ity c lass C^ {U) with 
respect to the zeroth-order derivative level of g. Note especially that the jump conditions (3.25) and (3.26) are 
precisely of such a nature that no violations of the jump conditions (3.14) and (3.15) above for the semi-longitudinal 
Weyl curvature characteristic eigenfields may occur. 



IV. WORKED EXAMPLE: COSMOLOGICAL MODELS WITH ABELIAN G2 ISOMETRY GROUP 



In this section we turn to discuss in some detail a new, fully gauge- fixed, (?^g-order autonomous FOSH evolution 
system for spatially inhomogeneous perfect fiuid cosmological models ( M, g, u ) which arc invariant under the trans- 
formations of an Abelian G2 isometry group that is simply transitive on spacelike 2-surfaces. Thus, all geometrically 
defined field variables vary in one spatial direction only. A systematic approach to this class of cosmological 
models was brought forward some time ago by Wainwright in Refs. [Q and [^ , wherein the generic case was given 
the classification label "A(i)". A number of exact, real analytic, solutions to the EFE for Abelian G2 perfect fiuid 
cosmological models are known, such as those listed in Refs. [^ and [^ or the singularity-free solution obtained by 
Senovilla [p5|, but most of them belong to dynamically restricted or higher-symmetry subcases. 



A. Well-posed Cauchy initial value problem 



Choosing an orbit-aligned group-invariant 1 + 3 ONF basis {bq, eo, } such that commutation relations 
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= [|,eo] = [^,ea] = [r7,eo] = [J7,e„] 



(4.1) 



hold between the two commuting spacelike Killing vector fields ^ and r) and { Gq, eo, }, and assuming that eo = u is 
orthogonal to the isometry group orbits, it follows that for all solutions in the Abelian G2 class we have ]30| , |3l|l 

= 62(7/^) = e3(w^) , = M2 = "3 '^"(u) = 02 = 03 = (n - 2n+) = n2 = . (4.2) 

That is, 62 and 63 are tangent to the isometry group orbits. Besides Gq = u also the frame basis field ei is hypersurface 
orthogonal |30|] . The canonical choice Wainwright proposes for Abelian G2 perfect fluid cosmological models consists 
of introducing fluid-comoving local coordinates {t, x, z} adapted to the integral curves of u, ^ and rj such that 

13 

^ = dy, v = d.; (4.3) 

the isometry group orbits, which have vanishing Gaufiian 2-curvature, are thus given by spacelike 2-surfaces {t = 
const , X — const}. Additionally, the coordinate components of the 1 + 3 ONF basis {eo, } as introduced by Eq. 



(2.10) are specialised to |31 



en = M ^dt , 62 = 62^ 9, 

ei = ei^ dx + ei'^ dy + ei^ dz , 63 = 63'^ dy + 63^ d, 



y ' - (4.4) 



in view of = w"(u) the choice = Mi is made.H Note that Wainwright's canonical choice establishes the property 
82 II ^. All geometrically defined field variables u will now only be functions of the local coordinates t and x. From 
the commutators (see, e.g., Refs. 19 1 and |l^), the canonical choice of 1 + 3 ONF basis and local coordinates has the 
direct consequences 

= (73 fix + r!i) = {V3(72 + ^2) = {V3a3 - ^3) , (4.5) 



implying the three spatial coordinate conditions = 62^ = e2^ = 63^ in Eq. (4.4) will be automatically preserved 
along u. Thus, the spatial frame { Gq, } will presently not be Fermi-transported along u, but fl" (and so the three 
frame gauge source functions T"/^) will be adapted to the fluid rate of shear instead. Likewise, the canonical choice 
leads to 

= (n+ - VSn^) . (4.6) 

We now define two new frame variables for two components of the fluid rate of expansion tensor by 

a:=(ie-2a+), /3:=(ie + a+), (4.7) 

shadowing the variable names used earlier in a related context in Ref. Q (on this choice of variables, see also the 
remarks made on the so-called "Taub gauge" for fluid spacetime geometries in Ref. |p8|). 

Substituting from Eq. ( [1.5| ) into the evolution equation for the semi-longitudinal fluid shear component CT3 (see, 
e.g., Ref. |[To|]), one flnds that the latter is involutive. Therefore, on performing a spatial rotation of { } about Gi 
at every point of a given 2-surfaces {t = const, x = const}, one can set 

= (73 (4.8) 

to hold at every event of (A^, g, u).^ Det ermining the evolution along u of the frame gauge source function 
{T^)~^T^i = ill as outlined in Subsec. p^I B| , at this point, from a dynamical viewpoint, the coordinate and spa- 
tial frame freedom has been completely fixed (modulo coordinate reparameterisation freedom given by f = i'(^)i 
= x'(a:), 2/' = y + /(x) and / = z + 5(2:) 111). 

It is now fairly straightforward to derive from the equations of a 1 + 3 orthonormal frame connection formulation 
of gravitational fields as given in Refs. p9| and pO| ] an evolution system of autonomous partial differential equations 
in FOSH format for the following set of eleven geometrically defined field variables = u^(t,xy. 



For = (j;"(u), the choice = Mi is the simplest one possible, but it is by no means compulsory. 
^Alternatively one could choose 63 || 77, and then set = (T2 instead. 
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(ei\ f3, (72, ai, fi, a, ui, (t_, tt-x, ctx, n_ )' 



(4.9) 



Note that the evolution of the frame coordmate components other than ei^ is decoupled from this set; ei^ itself forms 
a background field for pairwise dynamical interactions between {a, iii }, {ct-, Hx } and {cx, } in the sense of 
Geroch 



11- dimensional autonomous first-order symmetric hyperbolic evolution system: 







ivj f-'tei — 












Af-1 Of 3 = 


_ 3 fl2 _ 3 / 2 , 2 , 2 , 2 ^ , 3 2 1 /2 _ 

2/-' — ' X X/ 2 2 2V -L 




- A) 






M~-^ rt^rr^ — 
ivi Ut0 2 — 


\0 fJ y u — J U2 










M"^ dtoi = 


— B (ill + (2i ) H- 3 f/7- (J^ — 72^ (T \ 










M-^ dtfi = 


-(a + 2/3) (m + p) 






clM- 




cl e-i dxUi = 


[ - + ^2 _ 3 (-^2 _^ ^2 _ ^2^ _ ^2 ^) _ g^2 _^ ^2 


- a? - i (^i 4 


-P)] 


M- 


^ dtui - 


- cl e-i d^a = 


- a iti - (a + 2 /?) [ ^ + p) - ] ui 












— [2ai (a — /3) + 6 (n_ cix - nx (J-) ] 






M 




+ ei^ d^Uy, = 


- (a + 2 /?) cr_ + 2V3 CT^ - cri - 2^3 - {ux - 


2ai)nx 




M 




+ ei"^ dx(7- = 


— otrix — iixO- 






M 




~ ex dxU- = 


- (a + 2/? + 2\/3(T_)(Tx + (mi - 2ai - 2\/3nx)n- 






M 




- ex^ d^ax = 


— (a — 2-^/30--) + (ui + 2\/3nx) (Tx ■ 







(4.10) 
(4.11) 

(4.12) 
(4.13) 
(4.14) 
(4.15) 



(4.16) 
(4.17) 
(4.18) 
(4.19) 
(4.20) 

The characteristic condition for the set (|4.10| ) - ( 4.20 ) is invariantly given by = Q = {u°' (a)^ [ {— + 
clh^'')CbCc] [{—u'^u'^ + 1^ /i"^*^) Crf Ce where Ca :— Wa<P are the past-directed normals to characteristic 3-surfaces 
C:{4>{x^) = const}. The connection characteristic eigenfields associated with non-zero characteristic velocities are 
{ (ai CsMi)/(l + c2)i/2 } with w = ±cs and { (ct_ T rix)/V2, (ctx ± n_)/%/2} with v = ±1. 

Initial data which is invariant under the transformations of an Abclian G2 isometry group has to satisfy the following 
set of constraint equations: 

Initial value constraint equations: 

(4.21) 
(4.22) 
(4.23) 
(4.24) 

2 A (4.25) 
(4.26) 
(4.27) 
(4.28) 

the 3-Ricci curvature scalar *R of spacelike 3-surfaces S:{t = const} orthogonal to u is given by 

*R 2 (2 ei^ 9^ - 3 ax) ax-6{n'i+nl) . (4.29) 

The constraint equations ( 4.25 ) - ( 4.28 ) in this set are specialisations of Eqs. ( AlC ), (A5), (A8) and ( |A15| ), respectively. 
The remaining ones, Eqs. (4.21) - (4.24), derive from the commutators. 

Algorithm: Given an equation of state of the form p — p{fi), the initial data, which can be specified 
freely (modulo minimal differentiability requirements and the remaining coordinate reparameterisation 
freedom) as functions of the spatial coordinate a; on a spacelike 3-surface S:{t = const} orthogonal to 
the fluid 4-velocity field u, are the values of the variables {ex^, a, cr_, rix, (Tx, n_, /i}, together with 
the cosmological constant A. Specifying the values of the variabl es | 6 2^, £3^, 6 3"^, (3, a2, ax} at one point 
on this 3-surface, their spatial distribution follows from Eqs. ( |4.22 ) - ( 4.25 ), respectively, while ux is 



d,M = 
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= 


(C'com)l 


= ei"^ 93:62^ - (ai -f \/3nx) 


62^ 




= 


(C'com)2 


= ex^ dxC^i^ - [ax - \/inx) 




1- 2V3n- 62^ 
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(^£0111)3 
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(Cg) ■■= 
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[Cih ■■= 


ei^ dxP + ai (a — /3) + 3 (n 




- Hx 0--) 
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{Cxh 


(ei^ dx - "i ax + Viux) CF2 






= 


(Cpf)i : 


= ei^ + (/i +p) ui ; 
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determined through Eq. (4.28) and M through Eq. (4.21). Then all time derivatives of these variables 
are known. Note that one of the coordinate components ei^ and ei^ is arbitrary as functions of x, too, 
while the other follows from the unit-magnitude property of ei. It remains to specify boundary conditions 
for all variables to obtain unique solutions. 

Two subcases of importance are contained within the class of Abelian G2 perfect fluid cosmological models: the or- 



thogonally transitive subcase arises when 
subcase when additionally — 63^ ^ 

The initial value constraint equations ( 4.22 ) 
autonomous partial differential equations, where the characteristic speeds are | v 



ei = ei = (72, which itself specialises to the diagonal ("polarised") 
c Tx ^ n^, leading to a diagonal line element (see Refs. jl^ and [p2[). 

( 4.28 ) are propagated along u via a FOSH evolution system of 

0, according to 



7-dimensional autonomous constraint evolution system: 



dt{C,om)i =-(« + /?+ x/3 (7_) (C7com)l ~ e^" (Ci)i 

at(Ceom)2 = -{a + P- (Ccom)2 - 2V3ax (Ccom)i - (Ci)j 

M-lai(Ccom)3 = - 

M-^dt{CG) = - 
M~^dt{Ci)i = - 

dtiCpFh ^-2{a + l3) (Cpf)i ^ {a + 2 (3) [cl + c'^ ^ {fi + p)] (Cpf)i - 2 c^^ (^ + p) (Ci)i 



'(a + P) (CG)-4(ui+ai)(Ci)i 

- (a + 3/3) (Ci)i + (72 (Ci)3 -Kill- ai) (Cq) - i (Cpf)i 

- (a + 3/3 - V3a-) (^1)3 - 9^2 [Cih 

-2 d^V 



(4.30) 
(4.31) 
(4.32) 
(4.33) 
(4.34) 

(4.35) 
(4.36) 



The full set of propagation equations for the Abelia n G 2 perfect fluid cosmological models now forms a larger au- 
tonomous FOSH evolution system according to Eq. (2T), with the previous 11-dimensional system as a subset. The 
virtue of the larger system is that it explicitly shows that the set of dynamical field equations introduced is consistent: 
the initial value constraint equations are preserved by the time evolution equations (if they are true initially, they 
remain true thereafter). This completes the discussion on a well-posed initial value problem for this class of cosmo- 
logical models in the 1 + 3 orthonormal frame connection formulation of gravitational fields. 



We finally list the specialisations which the expressions ( A18| ) - (A23) for the Weyl curvature characterist ic cig en- 
fields undergo by the geometrical restrictions imposed by the Abelian G2 isometry group. We use Eq. ( 4.27 ) to 
eliminate derivatives dx<^2- 



Weyl curvature characteristic eigenfields: 



H- 

{E3 =F H2 
{E2 ± H3 

{E- T Hy, 

{By, ± 

Note that = (T2 ^ = (£'3 =F if2 



i ei^ d^ai - ^{a- (3) f3 - a'^ 



2n^ 



2ni 



1 ^2 
2 '^2 



- 2 (O-- - Tlx) (o-x + ?l_) + I (0-_ + Tlx) ((Tx - n_) 

- \/3 ((Tx ± n-) (72 

(/3 =F fli + VS (7_ =F V3nx ) (72 

T (ei^ dxT a - ai) (a_ + nx) ± 2\/3n_ ((7x ± n_) ± rix (/3 T ai) 
T (ei^ dxTa~ ai) {a^ ± 7i_) =F 2\/3 7i_ ((7_ ip nx) T J^- (/? + ai) 

{E2 ± H3) holds. 



(4.37) 
(4.38) 
(4.39) 
(4.40) 
(4.41) 
(4.42) 



B. Transport equations for jump discontinuities in outward first derivatives 

To give a graphic example, which, we believe, will also be of some interest in numerical investigations of dynamical 
features of Abelian G2 perfect fluid cosmological models, we conclude this section with a brief derivation of the 
transport equations that describe how physically relevant jump discontinuities in the outward first derivatives of the 
initial data are propagated along the bicharacteristic rays of the setting. We confine ourselves to modes with u 7^ 
relative to u. To this end, we first turn to Eq. (|]J) for each ofue{±Cs,±l}, leading to the conditions 

[d^ui] =±Cs [d^a] , [d^uy] = =F [d^cr^] , [d^n^] = ± [d^a-y] , (4-43) 
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respectively. This then gives from Eq. (2.5), together with the evolution subsystem (4.15) - (4.20), the relations: 



(i) V — ±Cs longitudinal modes: 

= (M-i dt T Cs ei^ 9, + i (M"^ dtCs T Cs ei^ d^c,) 



+ a+\{a + 2[5)[c-'^ -^{pL+p)~ cl]T\csUi±Csai)[d^{a±c-^ux)] ; (4.44) 



jump discontinuities in the outward first derivatives of initial data for {a, mi }, subject to Eq. (4.43), travel along the 
local sound cones. Note that, because of the general functional dependence Cg = Csip), the evolution of sound cone 
initial data typically leads to the formation of "shocks" . This phenomenon becomes impossible in the special case of 
a linear baryotropic equation of state with p{^) = (7 — 1) /i, 1 < 7 < 2, where = (7 — 1)^/^ ~ const. 

(ii) w = ± 1 transverse modes: 

= (M-i dt T ei^ 5, + a + /3 T Ml ± oi) [d^{a. T "x)] (4.45) 
Q^{M-^dtTei^d, +a + (5Tui±ai)[d^{ay.±n^)] . (4.46) 



jump discontinuities in the outward first derivatives of initial data for { (T_ , } and { CTx 7 } , subject to Eq. ( 4.43 ) , 
travel along the local light cones. 

V. CONCLUSION 

The main conclusion of this work is that in examining a set of dynamical equations for a physical system such as the 
relativistic gravitational field equations, completed by the equations for all needed auxiliary variables, the constraint 
equations are crucial in determining what information can be propagated along the characteristic 3-surfaces of the 
evolution equations when these are expressed in FOSH format. We have explicitly shown how to examine the constraint 
equations to determine whether jump discontinuities in the derivatives of the initial data can be propagated along the 
various characteristic 3-surfaces in the case of the relativistic gravitational field equations with a baryotropic perfect 
fiuid matter source, including the Weyl curvature variables. This makes it clear that such an investigation is needed to 
complement the determination of the set of characteristic 3-surfaces of any FOSH evolution system with any existing 
supplementary constraint equations, in order to determine which characteristic 3-surfaces in the set are physically 
relevant. 

The process outlined enables us to show why the characteristic 3-surfaces for the semi-longitudinal Weyl curvature 
characteristic eigenfields {E^ =F H2) and {E2 ± H-^) that are associated with \ v\ = ^, apparent in a straightforward 
reduction of order d^g of the evolution system of the relativistic gravitational field equations for a baryotropic perfect 
fluid matter source to FOSH format, cannot in fact be activated. It demonstrates that potentially associated semi- 
longitudinal gravitational radiation cannot occur, despite the occurrence of related characteristic 3-surfaces in the 
FOSH evolution system. It should be noted that the issue is not that the constraint equations are incompatible with 
the evolution equations, in the sense of not being conserved under the system's time evolution. On the contrary, the 
propagation of these constraint equations is indeed compatible with the existence of these modes, as can be shown 
by considering an extended FOSH evolution system that includes variables representing satisfaction of the constraint 
equations ||l^,|l^ (see also Refs. and |2^). The issue is that the constraint equations do not allow the setting 

of jump discontinuities in (the derivatives of) the initial data for [E^ =p H2) and {E2 ± H^) because the values of the 
components of the constraint equations themselves cannot suffer jump discontinuities: they have to be continuously 
zero from one spacetime event to any nearby one, i.e., everywhere. 

What this analysis does not do is to show in what manner the semi-longitudinal Weyl curvature characteristic 
eigenfields (£^3 =F H2) and {E2 ± -ffa), originally identified by Szekeres and shown there to have observable physical 



effects 26 1, will evolve in time, nor does it adequately characterise what freedom there is in setting initial data for 
these modes. It would be helpful to have some characterisation of the full freedom to assign these modes on an initial 
data 3-surface. 



The worked example for perfect fluid cosmological models with an Abelian G2 isometry group presented in Sec. IV 
features neatly the conceptual and mathematical advantages one can gain from combining the idea of FOSH evolution 
systems with a 1 + 3 orthonormal frame connection formulation of gravitational fields. It will be usable as a multi- facet 
test bed for numerical experiments of spacetime geometry evolution processes in relativistic cosmology and already 
provides the basis for work-in- progress on an interesting new scale-invariant, dimensionlcss dynamical formulation for 
the orthogonally transitive Abelian G2 perfect fluid cosmological models M . 
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APPENDIX: 



1. Suggestive example: linearised relativistic gravitational field equations in a metric approach with 

Regge— Wheeler coordinate gauge fixing 

Rooted in tradition, a metric approach to formulating relativistic gravitational dynamics is still more frequently 
encountered in the literature than, e.g., the (extended) 1 + 3 orthonormal frame formulation employed in this paper 
and in Rcf. This motivates the brief discussion of a metric approach based example derived from a paper by 
Kind, Ehlers and Schmidt ||l^ that highlights the issue of the physical role of the constraint equations in the selection 
of appropriate initial data sets. The example arises as a special subcase of linearised relativistic gravitational field 
equations that describe small adiabatic, non-radial perturbations of a star in hydrostatic equilibrium. Starting from 
a set of local coordinates {t, r, ip} and imposing Regge- Wheeler coordinate gauge fixing conditions, the Ansatz 
for the line element contains two real-valued metric functions / — f{t,r) and g = g{t,r)\ the angular behaviour of 
the line element can be given in terms of spherical harmonic functions Yim{'&, ^) (for more details see Sec. 2 of Ref. 
[p^). One can then obtain from the linearised Einstein field equations a second-order symmetric hyperbolic evolution 
system for / and g that takes the form 

- 9 + 9" = -;9' + ;|/ + 4j(l(i + l)-2)<,. (A2) 

Moreover, / and g are bound to satisfy the constraint equation 

= (^) ■■= a" - J /' + ' .9' - ^ (/ + 5) - J (/ - .9) • (A3) 



The characteristic 3-surfaces C determined by the principal parts of Eqs. ( |Al|) and ( |A2[ ) satisfy the conditions 



{t^r) ~ const. If we assume that Eqs. (Al) and (A2) hold for / and g, then we obtain for (C) the evolution equation 



= - (C) + (Cr + - {Cy ~ (C) . (A4) 

Hence, if on a spacelike 3-surface S: \t — const} real- valued initial data Iq^.^ ~ { /, f,g,g} satisfy = (C) ~ {C), 



then any solution of Eqs. (Al) and (A2) satisfies — (C) everywhere. 

It is instructive to note that from Eq. ( [A^ ) we can now add any arbitrary real- valued multiple a g" to Eq. ( |A2| ) to 
obtain a principal part of the form 

-g + {l + a)g" , 

which instead defines characteristic 3-surfaces C that satisfy the condition ( (1 -I- a)^/^ i + r ) = const. In particular, 
we can choose a = — 1, leading to r = const, which corresponds to a propagation speed \ v\ — 0. This implies that 
no information is transported from one event to a spatially separated nearby one. A different view considers the 
differentiability properties of / and g: if one chooses real-valued initial data such that { /, /', g, g' } are continuous^ 



then the constraint equation (|A3D im poses no condition on the differentiability of /" (and so on /, from Eq. (Al) ), 
while g" (and so g, from Eq. ( [A2[ ) ) is required to be continuous too. Hence, jump discontinuities are physically 
allowed only in the initial value of the former quantity and are propagated via Eq. (Kl[). 
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2. Constraint equations component- wise 



Fluid shear divergence equations/ Jacobi constraint equations: 

= -l (Ci)i = ei(i e + a+) - ^ (e2 - 3 aa + VSna) (ag) - ^ (eg - 3 ag - VSng) (as) 

- i (e2 + 2 U2 " a2 - \/3 n2) (wg) + i (eg + 2 Ug - ag + a/3 ng) (W2) 

— 3 ai C7-|_ + 3 (n_ ax — Ux <J_ ) + i (n — 2 n+) wi (A-5) 

= (Cj)i = ^ (ei - 2 ai) (n - 2n+) + (es - 2 as) (ag + VSna) + (eg - 2 ag) (as - V3n2) 

+ 2(ie-2(T+)a;i + 2(\/3o-2+f22)w3 + 2(\/3a3-03)a;2 (A6) 

= -^[{Ci)2 T (Cj)3] = ei(crg T - ^ wg T ^ 02) +62(0-- T^x ± ^ai) + eg(cTx ±n_ + -^i^i) 
+ ^ (62 - 3 a2 + 3\/3 712) ((7+) T (eg - 2 ag) (n + n+) - ^ 82(6) 

- 3ai (cTg =F "2) - 3 02 (cT- =F «x) - 3ag (ctx ± ) - (n+ - VSu-) 02 

— rix (\/3o-3 ± ^2) — n2 {V^a- ± ai) + -\/3ng CTx ± as J^- 



T ^ (2%/3 0-2 T 2 ug - 2 f}2 ± ag ± V3 ng) T ^ (2^3 CTx + 2f]i±in±in+±V3n_) 



W2 



T ;^(|e + 2CT+-2V3a_±2MiTai± V3nx)wg (A7) 

= [ (Ci)g ± (Cj)2 ] = ei((72 ± ng + ^ a;2 T 013) + e2((Tx ± n_ - ^ wi) - eg(<T_ =f "x T ai) 
± ^ (62 - 2 02) (n + n+) + ^ (eg - 3 ag - sVSns) (a+) - ^ eg(e) 
- 3 ai (1J2 ± ng) - 3 a2 (ctx ± n_) + 3 ag (cr_ Trix) + {n+ + \/3 n_) ag 
+ nx (\/3(T2 ± ag) - ng (\/3(T_ =f ai) - \/3n2 ctx ± 02 n_ 

± ^(2\/3(7gT2i(2 + 2f2g±a2T\/3n2)a;i± ^ (| 9 + 2cr+ + 2\/3c7- ± 2ui =f ai T \/3nx) ^2 

± ^(2\/3£7x -2niT^nT^n+±\/3n-)a;3 . (A8) 

Fitiid vorticity divergence equation: 

= (C2) = (ei - ui - 2 ai) (wi) + (ez - W2 - 2 02) (^2) + (eg - wg - 2 ag) (wg) . (A9) 
Generalised Gaufi-Friedmann equation: 



= (Cg) = 2 (2 ei - 3 ai) (ai ) + 2 (2 62 - 3 02) (02) + 2 (2 eg - 3 ag) (ag) 
1 2 _ 2 

6 " 3 



1I2 22 e/2i 2, 2, 2^ 

+ gn — |n+ — 6 (n_ + + 712 + ng) 



+ 6 (i e - a+) (i e + a+) - 6 (a2 + a2 + ^ + a^) 

+ 2 (wi - 2 f^i) wi + 2 (w2 - 2 ^^2) W2 + 2 (wg - 2 fig) wg - 2 - 2 A . (AlO) 
Weyl curvature divergence equations: 

= - i (C4)i = (ei - 3ai) {E+) - ^^ (62 - 3a2 + VSnz) (i;g) - ^ (eg - 3a3 - y/Sn^) (^2) 
+ iei(/^) +3n_ i;x - 3nx + 3(7- i/x - 3crx i?- 

- 3a;i if+ - I ((72 - V3 a;2) i?3 + § {<J3 + VStoa) H2 (All) 

= - i (C5)i = (61 - 3ai) {H+) - ^ (62 - 3a2 + \/3n2) (ifs) - ^ (eg - 3ag - \/3ng) (ifz) 

- i (/i+p)a;i +3n_ilx - 3nx H_ - 3cr_ iJx +3crx E_ 

+ 3 wi £;+ + I (c72 - \/3 a;2) Sg - | (ag + \/3a;g) £^2 (A12) 
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= 73 [ (<^4)2 T [C^h] = (ei ± 3a+ - 3ai) (^3 T ^^2) + (es T \/3a3 ± 3a;3 - Saa - VSn^) {E_ t H^) 
+ (63 T \/3(T2 T 3a;2 - Sas + \/3n3) {E^ ±H_) + ^ (ea T 3\/3(T3 t 3a;3 - 3a2 + 3\/3n2) (£;+) 
T ^ (es T 3\/3a2 ± 3a;2 - 3a3 - 3\/3n3) (//+) - 3^ e^in) T^{^Ji+p)^^ 

± \/3(a_ T nx) (£^3 ± H2) ± \/3 (ctx ± n_) (^2 T H^) T (3a;i ± n+) (£2 ± H^) (A13) 

0= ^[(C4)3±(C5)2] = (ei±3a+-3ai)(S2±if3) + (e2T V3(T3±3a;3-3a2- V3n2)(i;x ±i^-) 
- (eg T VScTa T 3^2 - Sag + VSng) (i^- T ^x) ± ^ (e2 T 3\/3cr3 T 3^3 - 3a2 + 3^3 712) 
+ :;75 (63 T 3^3 0-2 ±3 W2 - 3 03 - 3%/3n3) (£;+) - ^e3(/z) ± (/i + p)w2 

T \/3(a_ T nx) (^^2 T H^) ± (ax ± n_) (£;3 ± H2) ± (3a;i ± n+) {E3 t H2) . (A14) 
Momentum conservation equations: 

0= (Cpf)i =c^ei(/z) + (At + p)ui (A15) 
0= (Cpf)2 =c2e2(/x) + (Ai+p)u2 (A16) 
= (Cpf)3 = cl 63 (At) + {n + p)u3 . (A17) 

H^et/i curvature characteristic eigenfields: 

E+ = -l 61 (ai) + ^ (62 - 3\/3n2) (02) + 5 (^3 + 3\/3n3) (03) 

+ ^ (62 - a2 - ^ n2) (n2) - ^ (63 - 03 + 75 ^3) (ns) 
+ (i e + 0-+) cr+ + i (n - 2 n+) n+ 

- (cr_ - nx) (cr- + nx) - (o-x + ?^-) (c^x - n_) 
+ i (0-3 - ^2) (0-3 + "2) + 5 (0-2 + ns) ((72 - n3) 

+ n?. + rix + 5 ("^i ~ 2 ^1) '^i - 1 (1^2 - 2 ^2) W2 - g (^3 - 2 1^3) 0)3 

+ |(C?g)ii (A18) 

i?+ = - ^ (62 - 02 - \/3n2) ((T2) + ^ (63 - 03 + \/3n3) ((T3) 

+ i (61 + 2mi + fli) (wi) - i (62 + 2u2 + a2 - 3-^3 n2) (^2) 

- i (63 + 2 U3 + a3 + sVSns) (^3) 

- i (n- 2n+)(T+ - I (c7_ - nx) (ctx +n_) + | (a_ +nx) (ctx - ri-) 

- 5 (^3)11 (A19) 

{E3 ip iJ2) = T 1 61 (0-3 T n2 - ^ a;3 T ^ 012) ± i e2(a-_ T "x ± ^ ai) 

± i 63(fTx ± n_ + ^ wi) T ^ (62 - (i2 + \/3 n2) (a+) + ^ (63 - 2 03) (n+) 

- \/3 (ct2 ± ^ 03 T I "3) (ctx ± n_) - \/3 (0-3 ± ^ 02 ± 1 712) (cr_ T rix ) 
+ (i 9 + cr+) 0-3 ± I ai (0-3 =F 2 712) ± i (n - n+) 02 + | (n - 2 n+) 713 

T ^ n_ ((T2 T n3 ± ^ 03) ± ^ rix (0-3 ± n2 ± ^ 02) 
± ^ (2 i(3 =F W2 ± 2 Q2 + 03 + \/3 ri3) Wi 

- ^ (t^i - 2^1 =F ± %/3n_) i:i>2 ± ^ (2mi + ai - VSrix) W3 

- ^[(Cg)i2±(C3)3i] (A20) 

(£2 ± if3) = T 1 6i(£72 ± n3 + ^ a;2 T ^ 03) ± i e2(o-x ± n_ - ^ wi) 

T 5 63(ct_ T rix T ^ oi) - ^ (62 - 202) (n+) T ^ (63 - 03 - \/3n3) (ct+) 



17 



+ \/3 (cr2 ± ^ as T 1^3) (cT- T «x) - \/3 (0-3 ± ± 1 712) (ctx ± n^) 

+ (i 6 + cr+) 0-2 ± 5 fli (0-2 ± 2 ris) =F i (n - n+) 0-3 + i (n - 2 n+) ^2 
T ^ ?^- (era ± "-2 ± ^ 02) T ^ "-X (0-2 T "-3 ± 03) 
T ^ (2u2 ± W3 T 2f^3 + a2 - V3n2)uJi 



=F :77^ (2ui + ai + \/3nx) 



2^3 

1 

\/3 



ti^2 



2^3 



UJ3 



(Cg)31 T (C3)i2 _ 



(A21) 



{E^ =F ^^x) = Tei(cr- T "-x) ± 562(0-3 T "-2 + ^t^3 ± ^ 02) T ^ 63 ((72 ± ^3 - ^ lj2 ± ^ 03) 

+ ^(a2±^a3±2n3)((T2±n3)- ^((73±75a2T2n2)((T3Tn2) 
+ (i 6 - 2 cr+) cr_ ± 3 nx cr+ ± i (n + 2 n+) CTx + ^ (?^ + 4 n+) n_ 
± ai (cr_ q= 2 nx ) + 5 02 "2 + 5 03 "3 ± ^1 

± ^ (2m3 =F W2 ± 2fi2 + a3 - \/3n3) 0^2 ± ^ (2 U2 ± a;3 T 2 il3 + 02 + \/3 712)^^3 



2\/3 



_ (Cg)22 — (Cg)33 ± 2 ((73)23 ] 



(A22) 



(Sx ± i?-) = Tei(crx ± "--) ± 562(0-2 ± 713 - ^CJ2 ± ^ 03) ± 5 63(^3 =F "2 + ^Ci^3 ± ^7^02) 

- ^(a2±^a3±2n3)(a3Tn2)- ^(^T3±;75a2T2n2)((72±n3) 
+ (i 6 - 2 0+) o-x T 3 n_ 0-+ =F 5 ("- + 2 n+) o-_ + i (n + 4 n+) 
± ai (ctx ± 2 n_) — i 02 713 + i 03 7i2 ± Jt-x 

T ^ (2m2 ±a;3 =F 2fi3 + + \/3n2) W2 ± ^ (2 U3 =F ^2 ± 2 f72 + 03 - \/3n3) W3 



2V3 



[2 (Cg)23 T (C3)22 ± ((^3)33] 



(A23) 
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TABLE I. Conventions for (1 -I- 1 -I- 2)-decomposition 



Outward frame derivative (across J): 


ei 


Tangential frame derivatives (along J): 


62, 63 


(Semi-) Longitudinal tensor/connection frame components (wrt. ei): 


a+, a2, as / n, n+, n2, 713 


Transverse tensor/connection frame components (wrt. ei): 


a-. Ox 1 nx 
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